Runaway evaporation for optically dressed atoms 



David Wilkowski 

Institut Non Lineaire de Nice, Universite de Nice Sophia-Antipolis, CNRS, Valbonne 06560, France, and 
Centre for Quantum Technologies, National University of Singapore , Singapore 117543, Singapore. 

Forced evaporative cooling in a far-off-resonance optical dipole trap is proved to be an efficient 
method to produce fermionic- or bosonic-degenerated gases. However in most of the experiences, 
the reduction of the potential height occurs with a diminution of the collision elastic rate. Taking 
advantage of a long-living excited state, like in two-electron atoms, I propose a new scheme, based on 
an optical knife, where the forced evaporation can be driven independently of the trap confinement. 
^— ^ ' In this context, the runaway regime might be achieved leading to a substantial improvement of the 

, cooling efficiency. The comparison with the different methods for forced evaporation is discussed in 

the presence or not of three-body recombination losses. 

o, 

Quantum degenerate gases are now routinely produced in tens of laboratories across the world. The current meth- 
ods, even if they differ in their approaches, rely, on their final stage for forced evaporative cooling in a conservative 
trap. Two classes of traps; a high-field-seeking magnetic trap and a far-off-resonance optical dipole trap, are used 
independently, or in combination. In addition to design constraints, such as optical access, robustness and repro- 
i— T ducibility, the experimental apparatus should produce a degenerate gas with the largest number of atoms and the 
43 , best duty cycle. The new evaporation scheme, I describe in this paper, goes in that direction. 

Forced evaporative cooling in magnetic trap is usually performed using a RF probe or knife to coupled a trapping 
state to an anti-trapping one. The frequency of the RF knife is swept to lower the energy barrier of the truncated 
potential. During forced evaporation the confinement of the ultracold gas remains almost unchanged. In an optimized 
scenario, when the evaporation occurs on the whole surface where the RF resonance is fulfilled, the spatial density 
increases faster than the inverse of the mean velocity leading to a net increase of the elastic collision; the so-called 
c/3 ' runaway regime El- 

In magnetic traps, the atomic magnetic moment is frozen to fulfill the trapping condition. Spin degrees of freedom 
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can be released in far off-resonance optical dipole trap. In this case, the interplay between many-body systems and 
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atomic internal degrees of freedom open the door of a rich variety of new phenomena in ultracold atoms physic. In 
particular, this approach has been used to produce spinor condensates [2|-|j], fermionic spin mixtures and molecular 
ultracold samples [5l4l2j|. 

For some atomic species, evaporation toward quantum degeneracy is difficult in a magnetic trap. In those cases, 
evaporation in an optical dipole trap is an alternative. Elements, like Cs 0, 13 or 85 Rb EH, for example, have 
unfavorable collision properties and their scattering length needs to be tuned thanks to a magnetically induced 

Feshbach resonance [16| . Two-electron atoms, which carry any electronic ma gne tic momentum in their fundamental 

level, have also to be cooled with all optical methods. Recently, BEC of 40 Ca U3j_ff Sr El El, m Yb and 170 Yb 
[U, Fermi sea of 87 Sr [22] and 173 Yb [H and degenerate mixtures of 88 Sr- 87 Sr (H, 176 Yb- 174 Yb, 173 Yb- 174 Yb [H 
* \ have been reported. 

In dipole traps no energy selective coupling to an anti-trapping state, similar to the RF knife, has been implemented 
so far. Evaporation takes place because hottest atoms, with a mechanical energy higher than the potential height, 
t-H . are escaping the trap. Thus, the dipole trap might fulfil at least two major requirements. First, the potential height 
has to be in the order of few thermal energy, ksT such that the evaporation rate is significant. Second, the spatial 
• i-H , confinement and the compression of the atomic gas have to be such that elastic collisions can take place leading to fast 
thermalization. Different strategies have been used to improve the preparation of the cold sample before evaporation: 
For alkali atoms, Raman E3 or polarization gradient cooling (26| have been implemented to lower the temperature. 
For two-electron atoms, low temperatures are achieved with Doppler cooling on the intercombination lines E34HH3 - 
29| . Improvement of the spatial density in the dipole trap is also reported using adiabatic compression by dynamically 
changing the trap geometry, for example with an additional dimple trap E3 or a Zoom lens to displace one of the 
beam waist position [26| . The forced evaporation is performed by dynamically lowering the potential height. 

In the original experimental realizations of the BEC [30| and the Fermi gas [5[ in all optical devices, the gas was 
confined in a cross dipole trap to insure a good 3D confinement. The forced evaporation was carried out by lowering 
the power P of the trap light field. As a consequence the trapping frequencies, 



u oc VP, (1) 

and the trap confinement are also reduced. In this situation and in sharp contrast with RF knife, the elastic collision 
rate decreases during the forced evaporation and no runway occurs [3 if . Thus, it is crucial to counteract the reduction 
of efficiency of cooling during evaporation with a very good starting elastic collision rate. 
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New optical based trap schemes have been successfully implemented to limit and event almost suppress the reduction 
of trapping confinement during the forced evaporation. Those realizations have the same underlying idea which consists 
in decoup ling the potential height driving the evaporation, and the trap confinement controlling the thermalization 
rate. In [32| a tilt trap was used where spin polarized atoms of Cesium are held in a fixed dipole trap with a 
superimposed varying magnetic field gradient. The extra constant force pulls the atoms out of the trap. In all optical 
schemes, different versions of the dimple trap have been explored. In those schemes, the trap is made of at least 
two independent laser beams: a tightly confining one for 2D trapping and a wider one to close the trap in the third 
dimension. In an ideal scenario, only the power of wider beam is reduced [17ll20|. Clement and co-authors [33| report 
3D confinement and runaway evaporation in a stubble misaligned version of the dimple trap which works in a similar 
manner as the tilt trap of [33 |. 

In this paper, I propose a new method for runaway evaporation in an optical trap. In contrast with the methods 
implemented so far successfully, evaporation is not based on lowering the trapping lasers power but on a potential 
truncated with an optical knife. The paper is organized as follow: In section U I define the general requirements to 
implement the forced evaporation with an optical knife in a dipole trap. On the basis of a thermodynamical quasi 
equilibrated state, I derive rate equations which govern the evolution of some macroscopic quantities of the ultracold 
gas in the trap (section [TTJ) . Different evaporation strategies can be put in place according to the strength of the 
three-body recombination loss rate. Finally, I will draw the conclusion of this work in section IIIII 



I. DIPOLE TRAP WITH OPTICAL KNIFE 



I consider a two levels system in the optical domain where the excited (ground) state, labeled |e) {\g)), has a 
radiative lifetime long enough to disregard any spontaneous emission. Like in other evaporation schemes, far off- 
resonance lasers insure the 3D confinement of the cold gas in the \g) state. It is crucial however that |e) has to be an 
anti trapping state. The \g) — ¥ |e) transition is driven with a quasi resonant optical field at a detuning S with respect 
to the bare frequency difference of the two states. In the RWA approximation, the eigenenergies are 




E ± = - 1±< 1 4£— ^ , (2) 



where s = V e (r) — S + V g (r) and p — (V e (r) — S) ■ V g (r) are respectively the sum and the product of the ground and 
excited state energies dressed with the far off-resonance trapping lasers. Q is the Rabi frequency of the quasi-resonant 
field. This field is supposed to give any extra light shift contribution. In the adiabatic regime, an atom, initially in 
the trap, i.e. in the \g) state, is resonantly brought into the anti trapping state |e) at a distance where the avoid 
crossing occurs (V g (r) = — V e {r) + <5). Hence forced evaporation might be achieved by sweeping the detuning of the 
quasi-resonant field from zero toward a blue value. 

For a practical implementation, two-electron atoms are a straightforward choice since one of the long living states 
in the triplet spectrum can be used as an excited state. The 3 Po state is in particular interest because of the absence 
of an electronic magnetic momentum. Thus the forbidden ISo — > 3 Pq transition can be almost decoupled to stray 
magnetic fields. For the odd isotopes, the So — » 3 Pq transition, due to hyperfine mixing, becomes weakly allowed to 



single photon direct excitation. The bare linewidth of the transition is in the range of few mHz [34j which prevents 
any spontaneous emission during evaporation and provides enough coupling strength among the two states. For the 
even isotopes, the ISo ~ > 3 Pq transition is strongly forbidden. In that case, the 3 Pq state needed to be mixed to 
another short living states as depicted in (3f| H(| . 

Beside the issues of getting favorable conditions for evaporation, one has also to keep in mind that the light shift 
due to the off-resonance lasers should create an anti-trapping state for |e). Supposing that \g) is a trapping state, this 
condition automatically gets fulfilled for a two-level system where as in real atoms light shifts values can be engineered 
almost at will according to the presence of the other excited levels. If far infrared traps have to be excluded because 
they give the same negative sign of the electrical polarizability for both states at play, then it seems, however, that 
this new constraint can be overcome in most cases. For Yb, a popular near infrared laser is suitable (37j where 
as for Sr a dipole trap made with a red laser might be an alternative. For the latter case, the 5s 2 Sq — > 5s6p 1 Pi 
and 5s5p 3 Po 5s6s 3 <Si transitions respectively at 461nm and 679nm give the major contribution to the states 
light shifts (see figure [T^,). An example of the potential, dressed with a quasi resonant knife laser at a detuning of 
5 = 2ir ■ 340kHz, is depicted on figure [TJd. More values of the parameters used in this example are given in the caption 
figure [TJ I emphasize however that Q — 2tt ■ 40kHz yields an energy splitting of the two adiabatic potential curve in 
the microKelvin range at the avoid crossing. It is important to get a substantial energy splitting because evaporation 
only occurs if the escaping atom stays in the adiabatic state. Landau-Zener tunneling rate to the neighboring state 
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FIG. 1: (a) Simplified energy diagram and states involved in forced evaporation of 87 Sr with an optical knife in quasi resonance 
with the 1 So — > 3 Po transition. The coupling laser has a power of 50mW focus on a 50/im spot size. The Rabi frequency is 
Q, — 27r-40kHz. The dipole trap is produced with 660nm crossed laser beams focus on a 50/im spot size. The total power is 1W 
for a mean radial trap frequency of 270Hz. (b) shows the energy radial distribution, in waist unit, of the potential in energy 
recoil unit. The dash blue (red) curve corresponds to the 1 So ground ( 3 Po excited) state in interaction with the trapping laser 
only. The plain curves indicate the adiabatic potential dressed with the optical knife. In this example, the detuning of the 
knife is S = 2vr ■ 340kHz. 
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TABLE I: Important parameters used in the thermodynamical model of evaporative cooling. Their expressions are given for 
a 3D harmonic trap. The exact expressions, in the second column, are extracted from reference |l|]. The third column gives 
approximated expressions and last column shows their validity domain within an error smaller than 5%. 



can be evaluated following reference [38( and it turns out that it can be disregarded as far as fi.fi > k^T. Ml has to be 
interpreted as the energy splitting of the two adiabatic states at resonance and fcf,T as the mean energy of an escaping 
atom [l|. 



II. EVAPORATION IN A FIX OR VARYING TRAP 



I have shown that an ultra-narrow (or clock) transition may be used as an optical knife for evaporation in a Dipole 
trap. The initial high confinement is then kept almost fixed, in contrast with the original evaporation techniques in 
which the trap height is reduced. Let us now compare the two methods, i.e. when the trap confinement is varying 
or remains fixed. For that purpose, I will establish some relations among quantities that characterize the efficiency 
of the evaporation. The starting point would be to set a kinetic equation of the dynamics of the ultracold gas in the 
finite depth trap for both cases. As far as the elastic collision rate insures a good thermalization of the gas during the 
evaporation, it has been showed that the distribution function in energy of the gas is well represented by a Boltzmann 
distribution truncated at an energy e t corresponding to the energy threshold at which the evaporation occurs 
Furthermore, a simple and efficient cooling ramp consists in keeping the parameter 

rj = € t /kbT (3) 

large and constant. Here, T is the temperature of the gas. Hereafter these conditions are supposed to be always 
fulfilled. Following references [l| and [3l[, a thermodynamical model of evaporation can be resumed in a set of rate 
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TABLE II: Experimental data and numerical results for Bosons. The simulation stops when D 
experimental data are given for a similar condition, n.a. means: data not available. 
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TABLE III: Experimental data and numerical results for Fermions. The simulation stops when T /Tf 
experimental data are given for a similar condition, n.a. means: data not available. 
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equations, 



N = -[ lv {n,T) + ll (n,T)]N 



(4) 
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governing the temporal evolutions of the number of atoms N in the trap and their temperature. Here, "f v (n, T) is the 
evaporation rate related to the elastic collision rate i y(n,T) — Ariva(see table HJ. A is a dimensionless parameter. Its 
value depends on the type -fermionic or bosonic- of the atom and on the number n s of spin states involve in the s-wave 
scattering process. For spin polarized bosons, A — 1. For unpolarized Fermions A = n *~ . v and a are respectively 
the thermal velocity and the elastic scattering cross section. The expressions of a and c are also given in table Q] 
7/(n,T) stands for the atoms loss rate of the three-body recombination. Other loss processes such as one-body and 
two-body losses, are neglected [39j |. The extra heating rate of the three-body recombination corresponds to the second 
term on the right hand side in equation ([5])- It comes from the fact that the recombination rate scales like n 3 and 
thus occurs most likely in the high density region of the cloud where the potential energy is lower than its mean value 
40]. 

Initially, I set ji(n,T) — and three-body recombination losses will be discussed latter on. After simplification of 
equations ([5]) and ((4]), it turns out that N and T are linked though the simple relation: 



T <x N a . 



(6) 
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From equation ([6]), one derives other relations among quantities of interest such as the peak spatial density no, the 
clastic collision rate 7, the phase space density D, or if it concerns Fermions, the temperature in Fermi unit T/Tp. 
These quantities also depend on the confinement of the gas and of the trap truncated energy e t ■ Taking into account 
the gaussian shape of the dipole trap, expressions can be derived from [4l| . However when n 3> 1, the atoms are 
present at the bottom of the potential where it can be approximated by an harmonic potential characterized by its 
mean oscillation frequency u>. In the fixed trap configuration, u) remains constant during the evaporation. In contrast, 
in the varying one where the total power of the trapping laser is reduced and with respect to the relations (TTJ), ([3]) 
and (O, one has ui oc . Both trap configurations can be modeled introducing a parameter 9 such that 

wociV e 7 ; (7) 

where 9 = 0(= 1) means that evaporation occurs in a fix (varying) 3D trap. Finally, from relations ©and ([7]), one 
sets: 

n ociV 1 -^ 1 -^", (8) 

7 cxiV 1 -( 1 -l e ) tt , (9) 

D oc jV-i-!(2-»)« (10) 

— OC jV-* + *(2-»)<\ (11) 

Tf 

Since a > 1, these relations clearly show some differences among the two configurations (9 = 0,1). Firstly, the 
runaway regime, i.e. an increase in 7, could be reached only if 9 = 0. As a direct consequence, the cooling time is 
shortest for a fix trap. Secondly, and without any surprise, D (T/Tp) increases (decreases) for both configurations. 
However, the 9 = configuration leads to a better efficiency of the evaporation. 

One has to keep in mind however that relations (|6"l lll[) are expressed only in a lossless trap ( 7 ;(n, T) = 0) and for 
•q being large and constant. Those assumptions are not necessary fulfilled in some experiments. Thus, comparison 
between the model predictions and the documented evaporative schemes can be carried out through only with great 
care. Nevertheless an attempt is summarized in tables |H] & Mil where the raw data of some successful experiments 
toward degenerate gases of two-electron atoms are compared with the results of the model presented in this work. 
As was already mentioned and discussed in the original papers [l8[, (l~9l | and also confirmed here, the experimental 
evaporation ramps for 84 Sr reasonably stick to the 9 = 1 model. In contrast for 40 Ca, 87 Sr, 173 Yb and 174 Yb, 
either of the two evaporation models matches with the experimental data. Here, the trap is made with two intensity 
independent dipole lasers to which more complex experimental ramps are done. It is worth mentioning however that 
experimental data might be constrained between the two trap configuration model, where the 9 = one stands to be 
the most efficient in term of evaporation time t evap and remaining atoms Nt . 

Lets us now analyze how the three-body recombination losses affect the previous results. Indeed, reference [13] 
addresses three-body recombination as a major loss mechanism. It has been also reported that three-body recombi- 
nation limit the elastic collision rate in a trap configuration close to the 9 = configuration [33J . Firstly, I recall that 
the three-body recombination event scales with n 3 . Thus equation ((8]) and tables ILTl fc IIIII indicate that the fixed trap 
configuration may be more affected than the the varying one. The three-body recombination atoms loss rate is given 
by 

ll (n,T) = L 3 (n 2 ), (12) 

with 

L 3 = mC-a 4 . (13) 
m 

Here, n; stands for the mean number of atoms lost per collision event. < C < 70 is a dimensionless factor that 
might also vary with the scattering length a (see [4(| and references therein) . Both values of n; and C are of crucial 
importance for quantitative evaluation of the three-body recombination losses, but unfortunately they are not well 
documented for two-electron atoms. In [13], the authors report L 3 = 3 • 10~ 27 cm 6 /s for 40 Ca which can be converted 
to n/C = 16.5. Plugging this value into the model, the efficiency of the evaporation is shown in figured The blue (red) 
stars correspond to the 9 = (= 1) configuration. One notices that at low scattering length the three-body losses 
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FIG. 2: Remaining atoms at Df = 1 after a forced evaporation for different scattering length values ranging form 10ao to 
2000ao. Each point correspond to a simulation at a given scattering length. For all points n = 10, Ni = 10 6 , = 20/iK and 
rii — 1.3 ■ 10 14 cm -3 leading to an initial phase space density of Di = 10~ 2 . The blue (red) curves correspond to 8 = (= 1) 
configuration where as the stars (dots) refereed to mC = 16.5(210) (see text). The dash lines disregard the loss induced by the 
three-body recombination (niC = 0) for the 6 = (blue) and 8 = 1 (red) configuration. The inset shows, as function of the 
initial atoms number for mC = 16.5, the critical values a c of the scattering length for which the two trap configurations give 
the same result. 



have a limited impact and the 9 — configuration remains the best option (for more details about the comparisons 
see the figure caption). However above a characteristic value of the scattering length, a c ~ 385ao in that case, it 
seems more relevant to use the varying trap where the spatial density remains moderate. Similar conclusion can be 
drawn if one uses the most pessimistic value niC — 3 • 70. In that case, the inversion of the efficiency among the two 
trap configuration is observed at lower value: a c ~ 105ao. The values of a c are non universal and may be figured out 
for any experimental realization. For example, in the inset of figure [5] is shown the variation of a c as function of Ni. 



III. CONCLUSION 



If three-body recombination losses can be disregarded, the dynamic of the evaporation, at r/ constant, is ruled by 
a set simple set of nonlinear equations for N and T. The forced evaporation in two configurations of the dipole trap 
has been compared: the varying trap where the laser power of the dipole trap is ramped down leading to a reduction 
of the spatial confinement and the fix trap where the spacial confinement is kept constant. I show that the runaway 
regime, characterized by an increase of the elastic collision rate, is reached only in the fix trap configuration. As a 
consequence the efficiency of the evaporation in term of remaining atoms and cooling time is significantly improved. 
However, the spatial density increases during the evaporation and three-body recombination losses have also to be 
considered. In this context I show that there exists a characteristic scattering length above which the varying trap 
configuration, with lower spatial confinement, becomes the best choose. 

A practical implementation of the fixed trap can be done using the 3 Pq long living excited state in two-electron 
atoms as an anti-trapping state. Just like with RF knife for magnetic traps, the dipole trap can be truncated using 
an quasi resonant optical field on the — ► 3 Pq transition. The forced evaporation is obtained sweeping the laser 
detuning. 

The discussion in this paper had been carried out considering the lower branch of the adiabatic potential (in blue 
figure [TJd). If the atoms are now transferred to the upper branch (in red figure [T|b), they are confined in a quasi 2D 
bubble-like optically-dressed trap. Similar types of traps have been already proposed [42[ and accomplished [43| with 
RF-dressed state in magnetic traps. In those, the interplay between the Zeeman sub-levels structure and the vectorial 
nature of the RF field creates some holes through which evaporation is done Q. In a bubble-like optically-dressed 
trap, such holes do not exist at least in the two levels scheme considered here. However, forced evaporation can still 
take place sweeping the dressed laser frequency in the opposite direction to the one used for the forced evaporation 
in the lower branch of the adiabatic potential. 
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